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1. Introduction 

The strange quark content of the nucleon (A/^l^^lA/^) is a fundamental quantity on the nucleon 
structure. It represents the effect of strange quark on the nucleon mass Mj^ 

m,{N\s-s\N) 

where is the strange quark mass. This parameter is also relevant to direct experimental searches 
for the dark matter, as one of its candidates, namely neutralino, may interact with the nucleon most 
strongly through the strange quark content. Therefore, (A/^l^slA/^) is crucial to assess the sensitivity 
of the experiments [|l|]. 

We calculate (A'^|si'|A'^) from lattice QCD simulations in this study. The relevant nucleon three- 
point function is purely composed of a disconnected diagram, whose computational cost is pro- 
hibitively high with the conventional simulation methods. We overcome this difficulty by using the 
methods of the low-mode averaging [H Hi and the all-to-all propagator [|4p. 

Another advantage of our study is that chiral symmetry is exactly preserved by using the 
overlap quark action. We point out in this article that explicit symmetry breaking with conventional 
lattice fermions induces operator mixing leading to a possibly large contamination in (A/^liilA/^). 



2. Simulation setup 

We simulate Nf=2 + 1 QCD with the Iwasaki gauge and the overlap quark actions. The lattice 
spacing is determined as a = 0.112(1) fm using the Q. baryon mass as an input. We simulate two 
values of degenerate up and down quark masses mud = 0.035 and 0.050 on a xNt = 16^ x 48 
lattice. The strange quark mass is set to = 0.080 and 0.100, which are close to the physical 
mass fUs^phys = 0.08 1 fixed from Mk- We also push our simulations to lighter masses m„j = 0.015 
and 0.025 on a larger lattice 24^ x 48 but with a single value of =0.080. The four values of m„rf 
cover a range of Mj^, ~ 290 — 520 MeV with a condition MjiL > 4 satisfied. We have accumulated 
about 50 independent gauge configurations at each combination of and m^.. 

We calculate the nucleon three-point function 

Clpt{y,t,rcAtAts) = l^LI^-' \^{N{l^,t,rc + ^t)ss{^,t,src+^t,)N{y,t,src))] 

x,z 

-{ss{z,t,rc+^ts)) tr, [T{N{x,t,rc+ ^^t)N{y,tsrc))] } , (2.1) 

where A'^ = e"'"^ {ulCY5db)uc is the nucleon interpolating operator, and ss is the strange scalar oper- 
ator on the lattice. We denote the temporal separation between the nucleon source and sink (scalar 
operator) by At (Af^). This correlator is calculated with four choices of the source location, namely 
tsrc = 0, 12, 24 and 36 with y = 0, and two choices of the projection operator F = r± = (1 it 74)72, 
which correspond to the forward and backward propagating nucleons. We then take the average 

C3pt{At,Ats) = ^£{c3+(y,f.rc,Af,Af,)+C3^p-(y,A^,-f„.„A^f-Af,^f-Af,)} (2.2) 

to reduce its statistical fluctuation. We also calculate the nucleon two-point function CzptiAt) in a 
similar way. 
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We calculate the disconnected strange quark loop from ss and its vacuum expectation value 
{ss) using the all-to-all propagator. It is expected that low-lying modes of the Dirac operator D 
contribute dominantly to low-energy dynamics of QCD, and we calculate this contribution exactly 

{D-%Ux,y) = I ^v(')(x)v(')(3^)t (Dv(') = A»v(')), (2.3) 

where the number of eigenmodes Ne is 160 (240) on our smaller (larger) lattice. The contribution 
of the remaining high-modes is estimated stochastically. We generate a single complex Z2 noise 
vector rj (x) for each configuration, and divide it into Nd = 3x4xNt/2 vectors rj (x), which have 
non-zero elements only for a single combination of color and spinor indices on two consecutive 
time-slices. The high-mode contribution is then estimated as 

{D-%i,u{x,y) = £ V/("'WT]W(j)t, (2.4) 

where i/a^^) (x) is obtained by solving 

Di/a(^)(x) = (1 - i^w)r](^)(x), = Iv(')(x)v(')(3;)^ (2.5) 

We note that, in this report, Csp, on 24-^ x 48 is calculated without the high-mode contribution in the 
quark loop, our measurement of which is currently underway. In our study in Nf = 2 QCD [Q, |6|], 
we observed that the result for (A/^l^^lA/^) is well dominated by the low-modes and does not change 
significantly by ignoring the high-mode contribution. 

We can improve the statistical accuracy of the nucleon correlators by using LMA. Let us divide 
C2pt into two contributions Ciptjow and Czptjugh- C2ptjow is the two-point function constructed only 
from the low-mode part of the quark propagator {D^^)iow, and C2pt,high is the remaining contribu- 
tion. We can average over all possible source points (y,?.sTc) for C2ptjow The piece representing the 
nucleon propagation in C^pt can be calculated in a similar way. In this study, we average over all t 
for tgrc but 16 spatial locations for y at each t^rc- 

In our previous study in Nf = 2 QCD, we observed that the smearing of the nucleon source 
and sink operators N is crucial to obtain a clear signal of C^pt at reasonably small At. We therefore 
construct N from the quark field with the Gaussian smearing 

( / (0 3 
q.mrM=Y^^[l + —Hj j q{y,t), //x,y = I(V< + ^x,y+j), (2-6) 

where the parameters are set to ft) = 20 and N = 400. This smeared operator is not gauge invariant 
but we fix the gauge to the Coulomb gauge. We use ^.„„^ for both C^pt and C2pt- 

3. Strange quark content at simulated quark masses 

We extract the strange quark content {N\ss\N) on the lattice from the ratio of Cspt and C2pt 

R{At,At.) = -— ^ imN). (3.1) 
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Figure 1: Ratio 7?(Ar= 1 3, Af^) at ot„^ = 0.050 and = 0.080 as a function of Afj. The vertical lines show 
the locations of the nucleon operators. The noisy high-mode contribution to the quark loop is omitted in this 
plot. 



Figure 2: Fit result R{At) at m„^/ = 0.050 and nis = 0.080 as a function of At. The horizontal lines show the 
result of a constant fit in At and its error band. 



In order to identify a plateau of R{At,Ats), it is helpful to consider the same ratio but approximated 
by taking only the low-mode contribution in the scalar loop. Figure |l| shows an example of the 
approximated ratio as a function of At^. The significant change near At^ ~ and At is due to a 
contamination from excited states. We therefore fit R{At , Af to a constant form in At^ = [5, At — 5], 
which is well separated from the nucleon source and sink. Note that this fit and the following 
analysis on 16^ x 48 are carried out using R{At,Ats) without the approximation. While we use the 
approximated ratio on 24^ x 48 as mentioned in the previous section, we confirm in the analysis on 
16^ X 48 that the high-mode contribution to the quark loop has only small effect to (A/^l^sjA/^). 

The result of the constant fit in terms of At^ is denoted as R{At) and is plotted as a function of 
At in Fig. The stability of R{At) at At > 12 suggests that these data are well dominated by the 
ground state contribution. We extract (A'^l^^lA'^) by a constant fit to R{At) in this region. 



4. Chiral extrapolation to the physical point 

In Fig. ^ we plot the matrix element (A^|55|A^) as a function of niud- Our data show a mild 



4 



Nucleon strange quark content in 2 + l-flavor QCD 



K. Takeda 




Figure 3: Strange quark content (A^|ii|A^) as a function of niu^. Circles and squares show data at = 0.080 
and 0.100, whereas filled and open symbols are obtained on 16^^ x 48 and 24^ x 48, respectively. We also 
plot the linear and constant fits by solid and dashed lines. 



dependence on both m„rf and m,, which can be well described by a Unear form 

{N\ss\N) =CQ + ci^udmud + c\^sms (4.1) 

with ;^^/d.o.f. ~0.06. While we also attempt a chiral extrapolation based on heavy baryon chi- 
ral perturbation theory the chiral expansion turns out to have a poor convergence up to the 
next-to-next-to-leading order at our simulated quark masses. We therefore use the linear chiral 
extrapolation. The systematic error due to the choice of the fitting form is estimated by compar- 
ing with a constant fit, which also leads to a reasonable value of ;^^/d.o.f. ~0.44. We obtain 
(A^l^^lA'^) = 0.086(81)(107) at the physical point, where the first and second errors are statistical 
and systematic, respectively. 

This bare matrix element is converted to the renormalization invariant parameter 

/, = «^ = 0,0B(12)(16), (4.2) 

We note that this is in good agreement with our previous estimate in Nf = 2 QCD ^ as well as 
our indirect determinations through the Feynman-Hellmann theorem in Nf = 2 and Nf = 2+ I 
QCD \§. 



5. Renormalization of 55 



In (4^), fr^ is expressed with the bare quantities assuming that the operator m^ss is renormal- 
ization invariant as in the continuum limit. This is, however, not the case if chiral symmetry is 
explicitly broken by the lattice fermions of the choice. To illustrate this, let us consider the renor- 
malization of the scalar operator ss in the flavor SU (3) symmetric limit for simplicity. Using the 
flavor triplet quark field i//^, ss can be written as a linear combination of the flavor singlet and octet 
scalar operators 

{ss)P''y' = i { {xj/YY''^' - \/3 (v/A V)^''^"} , (5-1) 
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where is a Gell-Mann matrix, and we put the superscript "phys" to the renormalized operators 
to distinguish them from bare lattice operators. The singlet and octet operators are renormalized as 



(5.2) 
(5.3) 



in renormalization schemes which respect chiral symmetry. The operator (ss)''''^^ is thus expressed 
as 



phys 



1 



{ (Zo + 2Z8 ){ss) + (Zo -Z8){uu + dd)} . 



(5.4) 



This implies that ss can mix with the up and down quark operators unless Zo = Zg. 

The difference Z^—Zg comes from disconnected diagrams, such as that in Fig. ^, contributing 
only to Zo. These diagrams are in fact forbidden by chiral symmetry: the quark loop with the 
scalar operator ss = srSl + siSr has to flip the chirality, while such a flip does not occur with the 
quark-quark-gluon vertices in the diagrams. This implies Zo=Zg, which also holds at finite quark 
masses in mass independent renormalization schemes. The renormalization of ss thus reduces to a 
simple multiplicative renormalization 



[ss 



.phys 



Zsss, Zs=Zs = Zo, 



(5.5) 



and the operator rrigss is renormalization invariant provided that chiral symmetry is exactly pre- 
served as in our study. 



When chiral symmetry is explicitly broken by a lattice fermion formulation, (5.5) is modified 



as 



phys 



(Zo + 2Z8) (5^^)''" + (Zo - Zs)iau + dd)'"' + ^ + 



(5.6) 



The second term represents the mixing with the light quark contents due to Zq^Z^. This may lead 
to a large contamination in {ss)P^^y\ because the light quark contents induce a connected diagram, 
the magnitude of which is much larger than that of the disconnected one. Previous direct calcula- 
tions with the Wilson-type fermions [|l^, [ll], in fact obtained rather large values /y^ ~ 0.3 -0.5 



compared to our result ( |4.2[ ) probably due to this contamination. 

The explicit symmetry breaking also induces the mixing of the flavor-singlet operator with 



lower dimensional operators, such as the third term in ( |5.6| ), through the renormalization ( |5.2| ). This 
contribution must be subtracted as a part of the vacuum expectation value of ss. Due to the cubic 
divergence, this results in a large cancellation, with which the calculation is potentially noisier. 



6. Conclusion 



We calculate the nucleon strange quark content in Nf =2+ \ QCD. Chiral symmetry is exactly 
preserved by employing the overlap quarks in order to avoid large contaminations from the operator 
mixing. We also use the LMA and the all-to-all propagator to accurately calculate the relevant 
nucleon correlators. Our result of fx^ is in good agreement with our previous study in Nf = 2 QCD 
as well as our indirect calculations for Nf = 2 and 3 using the Feynman-Hellmann theorem: all our 
studies consistently favor small strange quark content /r^ ssO.02. 
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Figure 4: Disconnected diagram contributing the renormalization of the flavor singlet scalar operator (cross). 
In higher orders, more gluons propagate between the quark loop and propagator on the bottom. 

Numerical simulations are performed on Hitachi SRllOOO and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 09/10-09). This work is supported in part by the Grant-in-Aid of the 
Ministry of Education, Culture, Sports, Science and Technology, (No. 20340047, 21674002 and 
21684013) and by the Grant-in-Aid for Scientific Research on Innovative Areas, (No. 20105001, 
20105002, 20105003 and 20105005). 
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